Quenching Spin Decoherence in Diamond through Spin Bath Polarization 
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We experimentally demonstrate that the decoherence of a spin by a spin bath can be completely 
eliminated by fully polarizing the spin bath. We use electron paramagnetic resonance at 240 gi- 
gahertz and 8 Tesla to study the spin coherence time T2 of nitrogen-vacancy centers and nitrogen 
impurities in diamond from room temperature down to 1.3 K. A sharp increase of T2 is observed 
below the Zeeman energy (11.5 K). The data are well described by a suppression of the flip-flop 
induced spin bath fluctuations due to thermal spin polarization. T2 saturates at ~ 250 ^s below 2 
K, where the spin bath polarization is 99.4 %. 
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Overcoming spin decoherence is critical to spintron- 
ics and spin-based quantum information processing de- 
vices [H 0]. For spins in the solid state, a coupling to a 
fluctuating spin bath is a major source of the decoher- 
ence. Therefore, several recent theoretical and experi- 
mental efforts have aimed at suppressing spin bath fluc- 
tuations I, i, i, 1 0, i, i . One approach is to bring the 
spin bath into a well-known quantum state that exhibits 
little or no fluctuations [l^, [Tl| . A prime example is the 
case of a fully polarized spin bath. The spin bath fluctua- 
tions are fully eliminated when all spins are in the ground 
state. In quantum dots, nuclear spin bath polarizations 
of up to 60% have been achieved [12, [HI. However, a 
polarization above 90% is need to significantly increase 
the spin coherence time [l^ ]. Moreover, thermal polariza- 
tion of the nuclear spin bath is experimentally challeng- 
ing due to the small nuclear magnetic moment. Electron 
spin baths, however, may be fully polarized thermally at 
a few degrees of Kelvin under an applied magnetic field 
of 8 Tesla. 

Here we investigate the relationship between the spin 
coherence of Nitrogen- Vacancy (N-V) centers in diamond 
and the polarization of the surrounding spin bath consist- 
ing of Nitrogen (N) electron spins. N-V centers consist 
of a substitutional nitrogen atom adjoining to a vacancy 
in the diamond lattice. The N-V center, which has long 
spin coherence times at room temperature fl5, 16], is an 
excellent candidate for quantum information processing 
applications as well as conducting fundamental studies 
of interactions with nearby electronic spins [H, [13, [ill 
and nuclear spins [H, [Ig. In the case of type-lb di- 
amond, as studied here, the coupling to a bath of N 
electron spins is the main source of decoherence for an 



N-V center spin [iSl . |2l| . We have measured the spin co- 
herence time (T2) and spin-lattice relaxation time (Ti) 
in spin ensembles of N-V centers and single N impurity 
centers (PI centers) using pulsed electron paramagnetic 
resonance (EPR) spectroscopy at 240 GHz. By compar- 
ing the values of Ti and T2 at different temperatures, 
we verify that the mechanism determining T2 is different 
from that of Ti. Next, we investigate the temperature 
dependence of T2. 

At 240 GHz and 8.6 T where the Zeeman energy of 
the N centers corresponds to 11.5 K, the polarization of 
the N spin bath is almost complete (99.4 %) for temper- 
atures below 2 K as shown in Fig. 1(a). This extremely 
high polarization has a dramatic effect on the spin bath 
fluctuations, and thereby on the coherence of the N-V 
center spin. We find that T2 of the N-V center spin is 
nearly constant between room temperature and 20 K, 
but increases by almost 2 orders of magnitude below the 
Zeeman energy to a saturation value of ~ 250 /xs at 2 K. 
The data shows excellent agreement with a model based 
on spin flip-flop processes in the spin bath. The observed 
saturation value suggests that when the N spin bath is 
fully polarized, T2 is limited by the fluctuations in the 
^^C nuclear spin bath. 

We studied a single crystal of high-temperature high- 
pressure type-lb diamond, which is commercially avail- 
able from Sumitomo electric industries. The density of 
N impurities is 10^^ to 10^" cm~'^. The sample was ir- 
radiated with 1.7 MeV electrons with a dose of 5 x 10^^ 
cm"'^ and subsequently annealed at 900 °C for 2 hours 
to increase the N-V concentration [2^ . 

Electronic spin Hamiltonians for the N-V {Hmv) and 
N centers {Hm) are. 



'Electronic address: ' susumu@iqcd.ucsb.edu| 



Hnv - D[{Sl 

NV gNV , 



NV\2 



1 



A 



NV cNV 



rN 



(1) 



2 



J^N — /-is* • 9 



Bn 



ANs^.I^, 



(2) 



(a) 



where /is is the Bohr magneton and Bq is the magnetic 
field. S^^ and are the electronic spin operators for 
the N-V and N centers and 7^ is the nuclear si 
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g" is the slightly anisotropic g-tensor of the N center. 
D = 2.87 GHz is the zero-field splitting due to the axial 
crystal field [1^. Due to the tetrahedral symmetry of 
diamond lattice, there are four possible orientations of 
the defect principal axis of the ^'^N hyperfine coupling of 
and A^^ . In the present case, = 114 MHz for 
the (lll)-orientation and A'^ = 86 MHz for the other 
three orientations [24]. For the N-V center, A^^ = 2.2 
MHz for the (lll)-orientation [2^. The nuclear Zeeman 
energy and the hyperfine coupling between the N-V (N) 
center and ^^C and the nuclear Zeeman energy are not in- 
cluded here. The energy states of the N-V and N centers 
are shown in Fig. 1(b). 

The measurement was performed using a 240 GHz 
continuous wave (cw) and pulsed EPR spectrometer in 
the electron magnetic resonance program at the National 
High Magnetic Field Laboratory (NHMFL), Tallahassee 
FL. The setup is based on a superheterodyne quasiopti- 
cal bridge with a 40 mW solid state source. Details of 
the EPR setup are described elsewhere [H, [1^ . No op- 
tical excitation was applied throughout this paper, and 
no resonator was used for either cw or pulsed experi- 
ments. Fig. l(c)-(f) shows cw EPR spectra at room tem- 
perature where the magnetic field was applied along the 
(lll)-direction of the ~ 0.8 x 0.8 x 0.6 mm"^ single crystal 
diamond. The applied microwave power and field mod- 
ulation intensity were carefully tuned not to distort the 
EPR lineshape. Five EPR spectra in Fig. 1(c) corre- 
sponding to the N center are drastically stronger than 
the remaining signals which indicates that the number of 
N centers dominates the spin population in the sample. 
The N EPR peaks show the slightly anisotropic g-factor 
g^ which gives = 2.0024 and g^ = 2.0025 - 6 and 
is in agreement with the reported g-anisotropy of type- 
Ila diamond f27|. As shown in Fig. 1(d), we also ob- 
served the much smaller N-V resonances which shows a 
line for the (lll)-orientation in the right side and three 
lines for the other orientations in the left side. An over- 
lap of the three lines is lifted because the applied Bq 
field is slightly tilted from the (lll)-direction. Based on 
the EPR intensity ratio between N and N-V centers, the 
estimated density of the N-V centers in the studied sam- 
ple is approximately 10^^ to 10^* cm^'^. EPR lineshapes 
of the N {\ms = -1/2, mi = 1) ^ |1/2,1)) and N-V 
{\ms — —1) |0)) centers are shown in Fig. 1(e) and (f) 
respectively. The N center shows a single EPR line with 
a peak-to-peak width of 0.95 gauss. On the other hand, 
the N-V center shows a broader EPR line (the peak-to- 
peak width is 2.36 gauss) due to the hyperfine couphng 
between the N-V center and the ^''N nuclear spins. The 
estimated hyperfine constant is 2 MHz, in good agree- 
ment with a previous report [23j . 
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FIG. 1: (a)Spins of the N-V and N centers at room tem- 
perature and at 8.56 tesla and 2 K. At room temperature, 
where up and down spins are nearly equally populated, the 
N spin bath polarization is very small and therefore, the spin 
flip-flop rate is high. At 240 GHz and 2 K, the N spin bath 
polarization is 99.4 % and the spin flip-flop rate is nearly 
zero. (b)Energy states of the N-V and N centers. The energy 
levels are not scaled. The states are indexed by \ms,mi). 
Transitions indicated by solid lines are EPR peaks used to 
measure the spin relaxation times Ti and T2. (c)cw EPR 
spectrum at 240 GHz at room temperature when the mag- 
netic field Bo is applied along the (lll)-direction. No optical 
pump is applied. The strongest five EPR peaks around 8.57 
tesla are from N centers, (d) N-V EPR peaks. The intensity 
ratio between the left-most N and the right-most N-V is ~ 80 
which corresponds to 120:1 population ratio between N and 
N-V centers respectively. Other impurity centers were also 
observed (not indicated). (e)N centers EPR for the transition 
of \ms = -1/2,771/ = 1) ^ |l/2, 1). (f)N-V centers EPR for 
the transition of Ittis — —1) ^ |0). 
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FIG. 2: l/Ti for the N-V and N centers as a function of 
temperature. Solid lines are the best fit of the spin-orbit 
phonon-induced tunneling model written by Eq. [S] Inset of 
the graph shows Ti versus temperature in a linear scale. 



The temperature dependence of the spin relaxation 
times Ti and T2 was measured using pulsed EPR. An 
echo-detected inversion recovery sequence (tt — T — 7r/2 — 
T—TT—T—echo) is applied for Ti where a delay T is varied, 
while a Hahn echo sequence (7r/2 — r — tt — r — echo) is ap- 
plied for T2 where a delay t is varied : 2Si] . The area of the 
echo signal decays as a function of the delay time T and 
2t for Ti and T2 respectively and therefore can be used 
to determine the relaxation times. For the pulsed EPR 
measurement, we used the \ms = — l,m/ = 0} <-> |0, 0) 
transition for the N-V center and the \ms ~ ^1/2, mj ~ 
1} ^ |l/2, 1) transition for the N center (Fig. 1(b)). 

The Ti for both the N-V and N centers was mea- 
sured from room temperature to 40 K. Below 40 K where 
the Ti is longer than 10 seconds, an accurate measure- 
ment proved impractical as the drift of the superconduct- 
ing magnet {'^ 5 ppm/hour) becomes nontrivial on the 
timescale of the measurement. The Ti is obtained by fit- 
ting a decay exponential to the recovery rate of the echo 
area t/o ~ ae"^/^^ . As shown in the inset of Fig. 2, the Ti 
of both centers increases significantly as the temperature 
is reduced. For the N-V center, Ti changes from 7.7 ±0.4 
ms to 3.8 ± 0.5 s. For the N center, Ti increases from 
1.4 ± 0.01 ms to 8.3 ± 4.7 s. To evaluate the tempera- 
ture dependence of the N center, we applied a spin-orbit 
phonon-induced tunneling model which is independent of 
the strength of a magnetic field [1^. The temperature 
dependence is given by the following, 



— ^AT + BT^, 



(3) 



where A and B are parameters related to Jahn- Teller 
energy and electron-phonon interaction [29l |. From the 
fit, we found A = 8.0 x 10"^ ^^jj 5 = 3.5 x 10"!° which 
are in good agreement with the values in Ref. [1^, and 
confirm a largely field-independent Ti relaxation. The 
temperature dependence of the N-V center also shows 
similar behavior. The Ti relaxation mechanism for the 
N-V center is beyond the scope of this paper [soj . 
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FIG. 3: (a) Echo area of the N-V center as a function of delay 
2r measured at room temperature and T = 1.28±0.1 K. Solid 
line shows the best fit by the single exponential. (b)l/r2 for 
the N-V and N centers versus temperature. The scale of the 
main graph is log-log. Solid lines are the best fit using Eq. 
The arrow shows the Zeeman energy of 11.5 K. The inset 
shows T2 versus temperature in linear scale which shows a 
dramatic increase of T2 below the Zeeman energy. 



We also investigated the temperature dependence of 
the spin coherence time T2 for the N-V center using a 
Hahn echo sequence where the width of the pulses (typ- 
ically 500-700 ns) was tuned to maximize the echo size. 
Fig. 3(a) shows the decay of echo area at room temper- 
ature and at r = 1.28 ± 0.1 K. These decays, which 
are well fit by a single exponential e"^"^/"^^ as shown in 
Fig. 3(a), show no evidence of electron-spin echo enve- 
lope modulation (ESEEM) effects from the ^''N hyperfine 
coupling [28| . This is due to the relative long microwave 
pulses and the nuclear Zeeman splitting at 8.5 T which is 
much larger than the ^"^N hyperfine coupling of the N-V 
center. Between room temperature and 20 K, we ob- 
serve almost no temperature dependence with T2 ^ Ti , 
(e.g. the T2 = 6.7 ± 0.2 fj,s at room temperature and 
T2 = 8.3 ± 0.7 fj,s at 20 K). This verifies that the mech- 
anism which determines T2 is different from that of Ti. 
Below the Zeeman energy (11.5 K), T2 increases drasti- 
cally as shown in the inset of Fig. 3(b). By lowering the 
temperature further, T2 increases up to 250 fis at 1.7 
K and doesn't show noticeable increase below 1.7 K. 

At high magnetic field, where single spin flips are sup- 
pressed, the fluctuations in the bath are mainly caused 
by energy-conserving flip-flop processes of the N spins. 
The spin flip-flop rate in the bath is proportional to the 
number of pairs with opposite spin and thus it strongly 
depends on the spin bath polarization [31|. At 240 GHz 
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and 2 K, the N spin bath polarization is 99.4 % which 
almost eliminates the spin flip-flop process. This ex- 
periment therefore verifies that the dominant decoher- 
ence mechanism of the N-V center in type-lb diamond 
is the spin-flop process of the N spin bath. Using the 
partition function for the Zeeman term of the N spins, 
Z = Es{!_i/2e"'^^''^""^''^ where /3 = l/iksT) and kg 
is Boltzmann constant, the flip-flop rate is modelled by 
the following equation [3l| , 

Tjr = C'Aras=-l/2^ms = l/2 + ^res 
-12 

C 

= ^i + eTzJT^{l + e-Tz./T) + r-^: (4) 

where C is a temperature independent parameter, Tze 
is the temperature corresponding to Zccman energy and 
Tres is a residual relaxation rate. We fit the T2 data for 
the N-V center using the equation above. The fit was per- 
formed with the fixed Tres = 0.004 (/is^^) corresponding 
to 250 /iS. This model fit the data well as shown in the 
log scale plot of Fig. 3(b). Tze = 14.7 ± 0.4 K obtained 
from the fit is in reasonable agreement with the actual 
Zeeman energy of 11.5 K. The result thus confirms the 
decoherence mechanism of the N spin bath fluctuation. 

The observation of the saturation of T2 ~ 250 fis also 
indicates complete quenching of the N spin bath fluctu- 
ation and a second decoherence source in this system. 
From previous studies [l^, [l^ , the most probable second 
source is a coupling to the ^'^C nuclear spin bath. In fact, 
T2 ~ 250 /iS agrees with an estimated decoherence time 
of ^'^C spin bath fluctuations p^ . 

Finally we investigate temperature dependence of T2 
for the N center at 240 GHz. No temperature dependence 



of T2 was observed in a previous pulsed EPR study at 
9.6 GHz We measured the \ms = -1/2, mi = 1) ^ 
|1/2,1} transition shown in Fig. 1(b) which can excite 
only 1/12 of the N center population while it is assumed 
that all N spins in this transition arc on resonance [2^ . 
The temperature dependence of T2 therefore shows the 
relationship between 1/12 of the N center and 11/12 of 
the N spin bath fluctuation. Similar to the N-V center, 
we found slight change between room temperature and 
20 K, i.e. T2 = 5.455 ± 0.005 /is at room temperature 
and T2 = 5.83 ± 0.04 us at 20 K, and then a significant 
increase below the Zeeman energy. Eventually, T2 be- 
comes 80 ± 9 /iS at 2.5 K. As shown in Fig. 3(b), the 
temperature dependence of T2 is similar to that of the 
N-V center. These facts support strongly that the de- 
coherence mechanism of the N center is also the N spin 
bath fluctuation. 

In conclusion, we presented the temperature depen- 
dence of the spin relaxation times Ti and T2 of the N-V 
and N centers in diamond. The temperature dependence 
of T2 conflrms that the primary decoherence mechanism 
in type-lb diamond is the N spin bath fluctuation. We 
have demonstrated that we can strongly polarize the N 
spin bath and quench its decoherence at 8 T and 240 
GHz. We observed that T2 of the N-V center saturates 
~ 250 fis below 2 K which indicates a secondary deco- 
herence mechanism and is in good agreement with an 
estimated coherence time dominated by ^'^C nuclear spin 
fluctuations. 
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